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Abstract. In this paper we discuss in detail the nonlinear equations of the mean-field 
approximation and their connection to the exact many-body Schrodinger equation. 
Then we analyze the mean-field approach and the nonlinear dynamics of a trapped 
condensate of weakly-interacting bosons. 



1 Introduction 

In the past few years many authors, working in different fields, have shown great 
interest in the so-called quantum chaos or quantum chaology, i.e. the signature 
in the quantal systems of the chaotic properties of the corresponding (h — > 
0) semiclassical Hamiltonian [1-4]. Incidentally, as stressed by Berry [5], "the 
semiclassical limit h — > and the long time limit t — > oo are not interchangeable 
- the origin of the (h, t~ 1 ) plane is mightily singular". 

The subject is very wide but, for reasons of space, we focus our attention only 
on the connection between the mean-field approximation and the onset of chaos. 
For a quantum system with discrete spectrum, dynamical chaos is possible only 
as a transient with lifetime tn, the so-called Heisenberg time, which scales as 
the number of degrees of freedom. Because i# can be very long for a many-body 
system, we suggest that the transient chaotic dynamics of quantum states and 
the related observables can be experimentally measured. Moreover, when the 
mean-field theory is a good approximation of the exact many-body problem, 
one can use the nonlinear mean-field equations to estimate the transient chaotic 
behaviour of the many-body system. As a specific example, we consider the 
dynamics of a trapped weakly-interacting Bose-Einstein condensate. 



2 Variational principle and mean— field approximation 



Let us consider a A-body quantum system with Hamiltonian H. The exact time- 
dependent Schrodinger equation can be obtained by imposing the quantum last 
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action principle on the Dirac action 

s= J dt< m\^§- t - H\m > , (i) 

where ip is the many-body wavefunction of the system. Looking for stationary 
points of S with respect to variation of the conjugate wavefunction ip* gives 

ih^i> = m ■ (2) 

As is well known, it is usually impossible to obtain the exact solution of the 
many-body Schrodinger equation and some approximation must be used. 

In the mean-field approximation the total wavefunction is assumed to be 
composed of independent particles, i.e. it can be written as a product of single- 
particle wavefunctions <f>j. In the case of identical fermions, ip must be anti- 
symmetrized [6]. By looking for stationary action with respect to variation of a 
particular single-particle conjugate wavefunction </>* one finds a time-dependent 
Hartree-Fock equation for each <j>j : 

lh §t^ >= Hj ' (3) 

where h is a one-body operator. The main point is that, in general, the one- 
body operator h is nonlinear. Thus the Hartree-Fock equations are non-linear 
(intcgro-)differential equations. These equations can give rise, in some cases, to 
chaotic behaviour (dynamical chaos) of the mean-field wavefunction. 



3 Mean Field Approximation and Chaos 

In the mean-field approximation the mathematical origin of dynamical chaos re- 
sides in the nonlinearity of the Hartree-Fock equations. These equations provide 
an approximate description, the best independent-particle description, which 
describes, for a certain time interval, the very complicated evolution of the true 
many-body system. Two questions then arise: 

1) Does this chaotic behavior persist in time? 

2) What is the best physical situation to observe this kind of nonlinearity? 

To answer the first question, it should be stressed that quantum systems 
evolve according to a linear equation and this is an important feature which 
makes them different from classical systems. Since the Schrodinger equation is 
linear, so is any of its projections. Its time evolution follows the classical one, 
including chaotic behaviour, up to tu- After that, in contrast to the classical 
dynamics, we get localization (dynamical localization). The Liouville equation, 
on the other hand, is linear in classical and quantum mechanics. However, for 
bound systems, the quantum evolution operator has a purely discrete spectrum 
(therefore no long-term chaotic behaviour). By contrast, the classical evolution 
operator (Liouville operator) has a continuous spectrum (implying and allowing 
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chaos). This means that persistent chaotic behaviour in the evolution of the 
states and observables is not possible. Loosely speaking, chaotic behaviour is 
possible in quantum mechanics only as a transient with lifetime tn [7,8]. 

The Heisenberg time, or break time, can be estimated from the Heiscnbcrg 
indetermination principle and reads 



tH*±, (4) 



where AE is the mean energy level spacing and, according to the Thomas-Fermi 
rule, AE oc h N , where N is the number of degrees of freedom, i.e. the dimension 
of the configuration space. So, as h — > 0, the Heisenberg time diverges as 

t H ~ h 1 -" , (5) 



and it does so faster, the higher N is [9]. We observe that the limitation to 
persistent chaotic dynamics in quantum systems does not apply if the spectrum 
of the Hamiltonian operator H is continuous. 

Concerning the second question, it is useful to remember that, in the thermo- 
dynamic limit, i.e. when the number N of particles tends to infinity at constant 
density, the spectrum is, in general, continuous and true chaotic phenomena are 
not excluded [10]. 

We have seen that the Heisenberg time tn is very large for systems with many 
particles. This fact suggests that the transient chaotic dynamics of quantum 
states and observables can be experimentally observed in many-body quantum 
systems. Moreover, when the mean-field theory is a good approximation of the 
exact many-body problem, one can use the nonlinear mean-field equations to 
estimate the properties of the transient chaotic dynamics. 



4 Nonlinear dynamics of a Bose condensate 

In this section we discuss the mean-field approximation and the nonlinear dy- 
namics for a system of trapped weakly-interacting bosons in the same quantum 
state, i.e. a Bose-Einstein condensate [11]. In this case the Hartree-Fock equa- 
tions reduce to only one equation, the Gross-Pitaevskii equation, which describes 
the dynamics of the condensate [12]. Nowadays, this equation is intensively stud- 
ied because of the recent experimental achievement of Bose-Einstein condensa- 
tion for atomic gasses in magnetic traps at very low temperatures (about 10~ 7 
Kelvin) [13]. 

The Hamiltonian operator of a system of N identical bosons of mass m is 
given by 

N Ti 2 1 N 

A = E ( - ^ + M'i)) + 2 E V (ri,rj) , (6) 

i—l ij— 1 
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where Vo(r) is an external potential and V(r, r') is the interaction potential. In 
the mean-field approximation the totally symmetric many-particle wavefunction 
of the Bose-Einstein condensate reads 

V>(ri, ...,r N ,t) = 4>{v ll t)...ct){v Nl t) , (7) 

where 0(r,£) is the single particle wavefunction. By using the quantum varia- 
tional principle for the Dirac action we get the equation 

»ftj^(r,t)= [-|lv 2 + yo(r) + (iV-l) J d 3 r'T/(r,r')|0(r',i)| 2 ]0(r,t), (8) 

which is an integro-diffcrcntial nonlinear Schrodingcr equation. If the bosons 
are weakly interacting, it is possible to substitute the true interaction with a 
pseudo-potential V(r,r') = g5 3 (r — r'), where g = 4Trh 2 a s /m is the scattering 
amplitude and a s the scattering length. In this way we obtain the so-called 
Gross-Pitaevskii (GP) equation 

ift J^(r, t) = [ - ^V 2 + V (r) + g(N l)|0(r, t)| 2 ] </>(r, *) • (9) 

We now consider a triaxially asymmetric harmonic trapping potential of the 
form 

V{v)= l - m ul(\\ X i + \ly* + \lz*) , (10) 

where Xi (i = 1,2,3) are adimensional constants proportional to the spring 
constants of the potential along the three axes. 

It has been shown, using a hydrodynamical approach [14], that in the strong 
coupling limit the GP equation has exact solutions which satisfy a set of ordinary 
differential equations given by 

^ai + Xfai = 9 - , i = l,2,3. (11) 

These nonlinearly coupled ordinary differential equations describe the time evo- 
lution of the widths a of the condensate wavefunction ip along each direction 
1 . Here, we have defined the coupling constant g = {2/tt) 1 / 2 {N — l)a s /ao, pro- 
portional to the condensate number TV and the scattering length a s . Note that 
t = uj t and do = (li/mwo) 1 ' 2 is the harmonic oscillator length. 

1 The same equations can be obtained by minimizing the Dirac action S with a trial 
mean-field wavefunction ip(ri, ...,rjv,t) = (j>(ri, t)...0(rjv, t), where 

* (r ' t] = (^dmm) 1/4 i n 3 cxp { ~^Sm + imx ^ ' 

with {x\,X2,X3) = (x,y, z). Ui and /?» are the time-dependent variational parameters 
[15]. 
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The three differential equations correspond to the classical equations of mo- 
tion for a particle with coordinates cr^ and Hamiltonian 



For g ^ this Hamiltonian is nonintegrable and thus generic. As is well known, 
intcgrable systems are rather exceptional in the sense that they are typically 
isolated points in the functional space of the Hamiltonians and their measure is 
zero in this space. If we choose at random a system in nature, the probability 
that the system is nonintegrable is one [16]. 

The small oscillations and the nonlinear coupling of these modes have been 
studied by Dalfovo et al. for Ai = A2 = 1 and A3 = \/8 (axially symmetric 
trap) [14]. One of us (L.S.) has recently calculated the mode frequencies of the 
low energy excitations of the condensate in the case of the triaxially asymmetric 
potential [17]. These excitations correspond to the small oscillations of variables 
a's around the equilibrium point, corresponding to the minimum of the effective 
potential energy of H. The eigenfrcquencies u> for the collective motion, in units 
of u) , are found as the solutions of the equation 



cj 6 - 3 (A? + X\ + A§) co 4 + 8 (Xj\ 2 2 + XjXj + A2A3) uj 2 - 20X^X1X1 = . (13) 



Near the minima of the potential the trajectories in the phase-space are 
quasi-periodic. On the contrary, far from the minima, the effect of the nonlin- 
earity becomes important. As the KAM theorem [18] predicts, parts of phase 
space become filled with chaotic orbits, while in other parts the toroidal sur- 
faces of the integrable system are deformed but not destroyed. The study of 
this order-chaos transition for the Bose condensate in the triaxially asymmetric 
potential is currently under investigation by our group. 

5 Conclusions 

The main conclusion of this paper is that the use of mean-field approximation 
leads to nonlinear equations. As a consequence, in some cases, the behaviour of 
the wavefunctions may be chaotic. 

As a specific example, the Bose-Einstein condensation for weakly-interacting 
trapped bosons has been discussed in great detail. 
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